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We propose an improved composite-boson theory of quan- 
tum Hall ferromagnets, where the field operator describes 
solely the physical degrees of freedom representing the devia- 
tion from the ground state. In this scheme skyrmions appear 
merely as generic excitations confined in the lowest Landau 
level. We evaluate the excitation energy of one skyrmion. Our 
theoretical estimation accounts for the activation-energy data 
due to Schmeller et al. remarkably well. 

Introduction: Low-energy charged excitations are 
skyrmions EJ in quantum Hall (QH) ferromagnets when 
the Zeeman energy {g*psB) is small. Skyrmions were 
initially considered as solutions of the effective nonlin- 
ear sigma (NLcr) model and later studied also in 
a Hartree-Fock approximation Their existence has 
been confirmed experimentally [ I ^) . The skyrmion has 
an intrinsic scale k, which is a function of the gyromag- 
netic factor g* (k, — > as g* — * oo and k — ► oo as g* — ► 0). 

However, the NLcr-model approach suffers from prob- 
lems. Basically, the NLcr model describes correctly per- 
turbative spin fluctuations , but there is no reason a 
priori that it is capable to treat nonperturbative objects 
such as skyrmions. We mention explicitly three points at 
issue. First, it is an unsatisfactory theory since the rela- 
tion between the microscopic wave function and the clas- 
sical solution of skyrmion is totally unclear. Second, it is 
an inconsistent theory because the skyrmion and vortex 
excitations are treated entirely different objects though 
their wave functions become identical as k — > 0. Third, 
the predicted skyrmion energy IE 10 is too large com- 
pared with experimental data |^,^. In particular, the ac- 
tivation energy A has been predicted that A = y/n/SE^ 
as g* — > 0, with Eq — e 2 /s£b the Coulomb energy unit. 
On the contrary, the observed value is consistent with 
A ~ as revealed in recent experiments fl. 

In this paper we propose a new scheme of skyrmions 
based on an improved composite-boson (CB) theory [ p"H . 
In this scheme the field operator describes solely the 
physical degrees of freedom representing the deviation 
from the ground state, and the semiclassical property of 
excitations is determined directly by their microscopic 
wave functions. Charged excitations are skyrmions and 
vortices, which appear merely as generic topological ex- 
citations confined within the lowest Landau level (LLL). 
Their excitation modulates the electron density p(x) ac- 
cording to the soliton equation (T^]?]] , 



— V 2 \np{x)-p{x) 



po = vQ{x), 



(1) 



where v is the filling factor, po is the average density 
and Q(x) is the topological charge density determined 



by the wave function. The main part of the excitation 
energy is the Coulomb energy induced by this density 
modulation. The skyrmion excitation energy is shown 
to vanish as g* — ^ 0. Our formalism is free from the 
difficulties mentioned above. Throughout the paper we 
use the natural units H = c = 1. 

Vortices: Leaving the derivation of the improved CB 
theory later, we summarize first its idea applied to the 
single-component QH state. We denote the electron field 
by ij)(x) and its position by the complex coordinate nor- 
malized as z — (x + iy)/2£s- Any state |6) at v = 1/m 
(m odd) is represented by the wave function, 



6[x] = (Oj^(asi) • • • il>(x N )\&) = u[z]& hN [x], 



(2) 



where &ln[x] = U. r<s (z r - z s ) m exp[- J2r=i \ z r\ 2 ] is the 
Laughlin function, and lo[z\ = lo(zi, Z2, • ■ ■ , zn) is an ana- 
lytic function symmetric in all N variables. The mapping 
from the fermionic wave function &[x] to the bosonic 
function u>[z] defines a bosonization. We call the under- 
lying boson the dressed composite boson and denote its 
field operator by f{x). We derive that 



<S v [x\ = (0\(p(xi) ■ ■ ■ <^(a;Ar)|6) = u[z]. 



(3) 



The Laughlin state is represented by & v [x] — 1. A typ- 
ical vortex state is by S v [a;] = J\ r z r, implying that 
(ip(x)) = z in the semiclassical approximation. This is a 
highly nontrivial constraint leading to the soliton equa- 
tion (|l|), where Q(x) = Q v (x) is the vorticity density, 



Q V (x) = —Eijdidjlnitpix)) = S(x). 
Zm 

The vortex carries the electron number, 



AiV 



d 2 x[p(x) - po] = -v, 



(4) 



(5) 



as follows from the soliton equation. 

Skyrmions: We next summarize the idea of an im- 
proved CB theory applied to QH ferromagnets. We de- 
note the spin component by the index a(=t,J.). Any 
state at v = 1/m is represented by the wave function 
similar to (^). Let us explicitly consider the case when 
the spinor component is factorizable, 



&[x]=U 



w T (z r ) 

LU^(z r ) 



Sln [x] ■ 



(6) 



The ground state is given by 0J'(z) = 1 and uj^(z) = 0. 
From a set of two analytic functions u> a (z) we construct 



1 



the complex-pro jective (CP) field n(x) whose two com- 
ponents n a {x) are 



n a (x) = 



LU a ( Z ) 



(7) 



This represents the most general skyrmion configuration 
jij. The normalized spin field, or the nonlinear sigma 
field, is defined by s a (x) = n(x)' l T a n(x) with r a the 
Pauli matrices. The skyrmion configuration is indexed 
by the Pontryagin number JIJ , whose density is 



Q p (x) = ^-e abc e ij s a {x)d i s b {x)&>s c {x). (8) 

57T 



The skyrmion excitation modulates not only the spin 
density but also the electron density according to the 
soliton equation ([!]), where Q(x) = Q p (x) is the Pon- 
tryagin number density (Q). 

The simplest excitation is given by one skyrmion with 
scale k sitting at x = 0, whose wave function is f|. 



&\x] 



n 



k/2 



©ln \x] 



(9) 



In the limit k — > the skyrmion is reduced to the vortex 
created in the spin-polarized ground state. The scale 
k is to be fixed dynamically to minimize the excitation 
energy. The Pontryagin number density (||) is calculated 
for the simplest skyrmion (|^) as 



<3» . I W£ 

7r [r z + {k£ 



b 



(10) 



The skyrmion carries the electron number (||) , as follows 
from the soliton equation ([j]) with Q(x) = Q p (x). 

Bosonization: We now construct the improved CB 
theory explicitly. We start with the kinetic Hamiltonian 
for planar electrons in external magnetic field (0, 0, —B), 

Hli = ^M j d2x ^( x K P * - lP y)( P * + iPy)1>(x), (11) 



where P» 



-id-j + eAj is the covariant momentum with 



Aj = ^EjkXkB. Here, £i 2 = -£21 = 1 and en = e 2 2 = 0. 

We define the bare CB field by way of an operator 
phase transformation of the electron field 4>(x), <j>(x) = 
e - l& ( x )i[)(x). The phase field 0(a;) is chosen to at- 
tach m units of Dirac flux quanta 2ir/e to each elec- 
tron via the relation, EijdidjQ(x) = 2irmp(x), where 
p(x) = ^(x)ip(x) is the electron density. When m is 
odd, 4>(x) is proved to be a bosonic operator. The covari- 
ant momentum for the bare CB field is = Pf, + df.® = 
—idk — eekjdjA(x). The bare CB feels the effective mag- 
netic field B c ff = e _1 V 2 A(x) , which vanishes on the ho- 
mogeneous state (p(x)} — po realized at v = 1/m. The 
auxiliary field A(x) is solved as 



A(x) 



d y ln 



2£r 



q(v), 



(12) 



with g(y) = p(y) — p . The bare CB is the one used in 



literatures 13 



We proceed to define the dressed CB field <p(x), 



ip(x) 



ip(x) 



(13) 



by dressing the bare CB with a cloud of the effective 
magnetic field described by A(x). Substituting ( |l3| ) into 
(|ll|), the kinetic Hamiltonian is transformed into 

H K = f d 2 x^{x)(V x - iV y )(V x + iV y Mx), (14) 

where we have defined ^(x) = ^ (x)e 2 - A ^ x \ with which 
p{x) = ^(x)\p(x) — ip$(x)ip(x), and 



Pj = - ( £ jk + iS jk )d k A(x). 



(15) 



Analyzing the Lagrangian density we find that the canon- 
ical conjugate of <p(x) is not iip^x) but iip$(x). A type 
of the field operator ( |l3| ) was first considered by Read 
]l4j and revived recently by Rajaraman et al. [^5j. 

We suppress the kinetic energy by imposing the LLL 
condition, 

(v x + iv y Mx)\e) = -±-JL<p( x )\G) = 0. (i6) 

Solving this condition we find that the iV-body wave 
function © v [x] is an analytic function as in ( p[) . It is 
an easy exercise to derive the following relation JL5|, 

p*(a!i) ■ ■■<p t (x N )\0) = e LN [x]^( Xl ) ■ ■■^{x N )\0). 

Because of this relation the function u)[z] in the wave 
function (||) is given precisely by the formula (|^). 

When the wave function is factorizable, &[z] = 
Y[ Lo(z r ), the one-point function is given by (ip(x)) = 
uj(z). Using (13) we may set 



(17) 



because p(x) = 4>(x)^ 4>(x) . It is easy to see that the 
Cauchy-Rieman equation yields the soliton equation (Q) . 
Hence, the soliton equation may be viewed as a semiclas- 
sical LLL condition for topological excitations [jl2|| . 

One might question the hermiticity of the theory |L5"| , 
since the covariant momentum (|l5|) has an unusual ex- 
pression. It is related with the fact that the canonical 
conjugate of ^(x) is iip$(x) = iip^ (x)e 2A ^ x K It implies 
that the hermiticity is defined together with the measure 
e 2A(x) ^ g^jj a measure has arisen since the transforma- 
tion (|l^) is not unitary. The covariant momentum Vj is 
hermitian together with this measure. 

Quantum Hall Ferromagnet: We proceed to ana- 
lyze the QH system with the spin degree of freedom. The 
electron field, the bare and dressed CB fields are denoted 
by ip a (x) , cj) a (x) and (p a (x) , respectively. We also use the 
two-component electron field ^f(x) and dressed CB field 



2 



&(x). All formulas from (11) to ( |l7| ) hold with appropri- 
ate modifications. The phase field Q(x) and the auxiliary 
field A(x) are defined by the same equations as before 
together with p(x) = *t(as)*(a;) = & (x)e 2A ^$(x). 
Now, the bare CB field is decomposed into the U(l) field 
4>(x) and the SU(2) field n a (x), (j) a {x) = <j)(x)n a (x) with 
n a the CP field, 



cp a (x) 



(x)n a (x). (18) 



The spin densities are 5* a (a;) = ^p(x)s a (x). The Zceman 
term is Hz — —g*psB J d 2 xS z (x). The ground state 



So 



is unique, upon which — 1 and 



,1 



0. Each 



Landau level contains two energy levels with the one- 
particle gap energy g* pbB. 

When the wave function is factorized as in (||), the 
one-point function is given by (ip a {x)) = u a (z). Based 
on the formula (fig) it is parametrized as 



(19) 



e- A ^e ix ^^/p~(xjn a (x)=u) a (z) ) 



from which the soliton equation ([!]) follows as the Caucy- 
Rieman equation 

On the other hand, the effective Hamiltonian describ- 
ing the perturbative property of the spin texture is de- 
rived as follows @,||. We make the LLL projection of 

the spin texture, which results in |S) = e l °\go) where 
O is a LLL projected SU(2) generator. We evaluate the 
energy of the state |6) and identify it with the effective 
Hamiltonian, H e g = (6\(Hc + Hz)\&), where He is the 
Coulomb interaction term. Making a perturbative ex- 
pansion in terms of the sigma field, making a gradient 
expansion and taking only the lowest order terms, we 
obtain [0,§ 

H cS = ^ Ps J d 2 x^d k s(x)} 2 - ^g*(i B Bs z (x)j, (20) 

as yields the NLcr model The first term represents 
the spin stiffness @ with p s = i/e 2 /(16\/27r££s). 

We consider the vanishing limit of the Zeeman term 
(g* = 0). In this case the ground state is given by an 
arbitrary constant sigma field, s(a;)=constant. All spins 
are polarized into one arbitrary direction. There exists a 
degeneracy in the ground states. The choice of a ground 
state implies a spontaneous magnetization, or a QH fer- 
romagnetism. When a continuous symmetry is sponta- 
neously broken, there arises a gapless mode known as the 
Goldstone mode. Quantum coherence develops sponta- 
neously. Actually, due to the Zeeman term all spins are 
polarized into the z axis. As far as the Zeeman effect 
is small enough, the system is still considered as a QH 
ferromagnet with a finite coherent length. The Gold- 
stone mode describes small fluctuations of the CP field 
around the ground state \go). We may set n" (x) = 1 and 
nHx) = C(aO/VPo, with iCHx),CM = i5(x - y). The 
effective Hamiltonian (p0|) implies |t]] 



H 



c -?£i i d 2 x[d k (Hx)d k ax)+& 2 (Hx)ax)}, (21) 



where £j, = y // 2p s /g*psBpo is the coherent length. 

Activation Energy: We analyze the soliton equa- 
tion (|l|) for the skyrmion excitation with the topological 
density (|l0|). Approximate solutions are constructed in 
the two limits, the large skyrmion limit (n 3> 1) and the 
small skyrmion limit (k <C 1). First, in the large limit 
we can solve (Xm iteratively, where the first order term is 



g(x) = p{x) - po 



\x), 



(22) 



with the Pontryagin density (|l0| ). This agrees with the 
formula due to Sondhi et al. j0f However, in the small 
limit the topological charge Q (x) is localized within the 
core: Indeed, we have Q p (x) — > 5(x) as k — > in (ITT 
with which the skyrmion is reduced to the vortex. 

We evaluate the excitation energy of one skyrmion. In 
the semiclassical approximation it consists of the electro- 
static energy Eq and the Zceman energy Ez, 



^skyrmion — Ec + -g* PbBAN s , 



where 



r e ' I ,2 ,2 Q(x)g(y) 

Ec = Ye I y ~ — ~- 



and the skyrmion spin AN S is 



AA, 



d 2 x{p Q - p(x)s z {x)}. 



(23) 



(24) 



(25) 



In evaluating the Coulomb energy we use ( |22[ ) for a large 
skyrmion to obtain Ec — v 2 ([3/k)E c with [3 = 37r 2 /64. 
For a small skyrmion we solve numerically the soliton 
equation (|l|) with Qix) = 5(x), and obtain Eq = v 2 ^E n c 
with 7 ~ 0.39. 

The skyrmion spin AN S would diverge logarithmically 
for the above skyrmion This is a fake since the the 
Zeeman term breaks the spin SU(2) symmetry explicitly 
and introduces a coherent length £l into the SU(2) com- 
ponent. The skyrmion configuration is valid only within 
the coherent domain because the coherent behavior of 
the spin texture is lost outside it. By cutting the upper 
limit of the integration at r ~ £l in (|2^) , we obtain 



AN, = k 2 In 



4t 

K i B 



1 



(26) 



with the coherent length £^ given below (|21|). 

We consider the integer QH state dXv — 1 . The energy 
of the skyrmion is explicitly given by 



^skyrmion - \ „ + ^ \8gK 2 



1 \&C, 



(27) 



where g = g* psB / Eq is the Zeeman energy in unit of 
the Coulomb energy Eq. By minimizing this energy the 
skyrmion scale k is determined as 



3 



/3 1/3 ^ln 



32<? 



-1/3 



(28) 



Our main result is the energy formula fl27|) with (28). It 
is notable that the skyrmion energy vanishes as ~g — > 0, 
where the skyrmion scale is infinitely large. 

The activation energy of a skyrmion-antiskyrmion pair 
will be given by 2i? s k yrm i on , which should be compared 
with experimental data. Our theoretical estimation is 
about two times bigger than the observed data for sample 
SI1 due to Schmeller et al. ||. There are two main rea- 
sons for the discrepancy: First, we have neglected a finite 
thickness of the layer. Second, we have assumed that the 
Coulomb energy of the antiskyrmion is the same as the 
skyrmion. Hence, the value f3 representing the strength of 
the Coulomb energy is subject to a modification in (27). 



We have fitted the date by adjusting the parameter (3 in 
Fig.[l]. Our theoretical curve reproduces the data remark- 
ably well. There is a small deviation for lj > 0.03, where 
the large-skyrmion approximation is no longer valid. The 
skyrmion spin is estimated that AN S ~ 8.7 at g = 0.01 
and AN S ~ 5.6 at g = 0.02. 

Discussions: We comment on the spin-stiffness term 
in (^0|) since our formula ( p3| ) is different from the stan- 
dard one |2|,[n| by this term. The term has been derived 
perturbatively from the microscopic Hamiltonian 
It describes correctly perturbative long-distance physics. 
However, there is no reason that it describes a nonpertur- 
bative object such as skyrmions. As we have argued, its 
absence is required from the consistency condition that 
the skyrmion wave function is reduced to the vortex wave 
function in the limit K — > 0. Indeed, the skyrmion energy 
is reduced to the vortex energy without this term. We 
also emphasize that the experimental data are excellently 
fitted by the skyrmion energy (|H) without it. A decisive 
test is the vanishing of the activation energy at g = in 
pure samples, which seems to be confirmed in a recent 
experiment |6). 
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